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Q , Abstract 

We consider the Cauchy problem to the two-dimensional incompressible liquid crystal 
I/"") ■ equation and the heat flows of harmonic maps equation. Under a natural geometric angle 

condition, we give a new proof of the global well-posedness of smooth solutions for a class 
of large initial data in energy space. This result was originally obtained by Ding-Lin in [3] 
and Lin-Lin- Wang in |11) . Our main technical tool is a rigidity theorem which gives the 
coercivity of the harmonic energy under certain angle condition. Our proof is based on a 
, frequency localization argument combined with the concentration-compactness approach 

which can be of independent interest. 
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1 Introduction 

m 

We consider the following hydrodynamic system modeling the flow of nematic liquid crystal 
\Q ■ materials in two dimensions (see, for instance, 0[71[8]): 

(N 

' u t + {u- V)u + Vp = An- V- (Vrf®Vd), 



in 
o 

(N 



d t + (u-V)d = Ad+\Vd\ 2 d, (1.1) 
V-n = 0, |d| = l, 



where u is the velocity field, p is the scalar pressure and d = (^1,^2,^3) is the unit- vector 
(\d\ = \J d\ + d| + d\ = 1) on the sphere § 2 C M 3 representing the macroscopic molecular 
orientation of the liquid crystal materials. Here the i th component of V • (Vg? <S> Vd) is given 
by Vj(Vjd • Vjd) = Y^j=i^jfflid ' Vjd). We apply Einstein's summation convention over 
repeated indices throughout this paper. For simplicity, we have set all the positive constants 
in the system to be one. We are interested in the Cauchy problem of (ll.ip with the initial 
data 

u(0,x) = uo(x), d(0, x) = do(x), lim do(x) = n G S 2 , (1.2) 

\x\— too 
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where n is a given unit vector. 

The above system (jl.ip is a simplified version of the Ericksen-Leslie model for the hydro- 
dynamics of nematic liquid crystals [5, 71 . The mathematical analysis of the liquid crystal 
flows was initiated by Lin and Liu in [51 [9]. In [15], Wang proved the global existence of 
strong solutions for rough initial data with sufficiently small BMO^ 1 norm (for uq and Vc?o). 
See also [10] for a small data global existence result in 3D. For any bounded smooth domain 
in M 2 , Lin, Lin and Wang [11] obtained the existence of global weak solutions which are 
smooth everywhere except on finitely many time slices (see also [6]). The uniqueness of weak 
solutions in two dimensions was studied by [121 [T6] . Very recently, a family of non-trivial 
two-dimensional exact large solutions was constructed in [3]. 

In this paper, we are concerned with the global existence of large classical solutions to 
(jl.ip . As always with these types of problems to extend the smooth local solution globally 
in time one needs to obtain certain a priori estimates. By the regularity theory in [ll|*l the 
smooth local solution (u, d) to (II. ip can be continued past any time T > provided that we 
have 

^ T (ll^,-)lll4 + ||Vd(t,-)Hl4)dt<oo. (1.3) 
The basic energy inequality associated with (jl.ip is the following: 

IflKt, -)|& + ||vd(t, OHIO + jT(l|Vti(., -)\\ 2 L , + ||Ad( S) •) + \vd\ 2 d( s , (i.4) 

< ^(Kllia + llVdoll^.), Vt>0. 

By (jl.3p ~ (jl.4|) and the Gagliardo-Nirenberg inequality 

i i 

||n|| L 4 (R2) < C ■ \\u\\ 2 L2{M2) ■ ||V«||£ 2(]K2) , 

it is obvious that 

i i 

H M HiL([0,T)xR2) < C ■ II^II2oo L 2 ([0iT)xR2) • ll^ll^flo.TjxR 8 ) 

< C • (|K|| L 2 + \\Vd a \\ L 2) < oo, for any T > 0. 
Hence the non-blowup criteria (ll.3p can be sharpened to 

/ ||Vd(t,-)|li4<ft < oo. (1.5) 
J o 

Note that in (jl.ip if we take u = we obtain the usual heat flows of harmonic maps. 
The condition (jl.5p asserts that the regularity of the whole liquid crystal system (jl.ip is in 
some sense governed by the pure heat harmonic flow part. For smooth solutions, the main 
obstruction in obtaining the a priori estimate (jl.5p is that the harmonic energy Eh arm {t) = 
||Ad(£) + |Vd(t)| 2 d(t)||2 in (jl.4p is not coercive. In particular it yields no useful lower bounds 
for general maps d : M 2 — > § 2 (a trivial example is the constant map). A natural idea is 

1 Although the main results in [TT] are stated for the bounded domains in E 2 , it is not difficult to check 
that the arguments there carry over to the R case. 
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to introduce some geometric constraints and work with a set of " restricted maps" such that 
the resulting harmonic energy is coercive. This idea was first used in the seminal work of 
Ding-Lin in [3] when they were studying a generalization of the Eells- Sampson's theorem. 
Inspired by [3], we follow this line of thought and establish global wellposedness for a family 
of initial data under a geometric angle condition. Roughly speaking, this condition says that 
the image of the initial orientation vector do is strictly contained in a hemisphere. The set 
of such maps is invariant under the dynamics of (jl.ip thanks to the maximum principle. 

Theorem 1.1. Denote e 3 = (0,0,1) G § 2 . Let u ,Vd G L 2 (R 2 ) with \d \ = 1, d - e 3 G 
L 2 (IR 2 ) and satisfy the geometric angle condition: 

inf d 03 > 0, (1.6) 

xeR 2 

where c?03 the third component of do. Then there exists a unique global smooth solution 
(d,u) to the incompressible liquid crystal equations (jl.ip with the initial data (|1.2p . Moreover, 
one has 

y°°(||Ad(i,-)||i 2 + ||Vd(t,.)||l4 + \\Vd(t,-)\\ 2 L2 )dt < oo. (1.7) 

Remark 1.2. We emphasize that the result in Theorem 11.11 (and also Theorem 11.71 below) 
can be obtained by using the method of [lljl 2 ! Our proof here is based on a frequency 
localization argument combined with the concentration-compactness approach which can be 
of independent interest. 

Remark 1.3. In Theorem 11.11 the choice of e 3 is for convenience only. In general one 
can choose any reference vector n G § 2 and require that the image of d is contained in a 
hemisphere around n. We should also point it out that the geometric angle condition (jl.6p 
may be necessary due to some counterexamples in heat flows of harmonic maps. In [1] the 
authors proved the global well-posedness of large solutions to the heat flows of harmonic 
maps for a class of initial data with symmetry: 

(xir^ 1 sin ^o( r )\ 
X2r~ l sin^o( r ) 
cosV'oW / 

where 

V>o(0) = 0, ip (R) < vr for all R > 0. 

Finite time singularities are also shown to exist in the case that vpo(R) > vr for some R > 
(see [2] for more details). 

Remark 1.4. As was shown in the previous remark, there exists a family of blowup solutions 
for the heat flows of harmonic maps when the geometric angle condition is violated. These 
blowup examples furnish a trivial (u = 0) set of counterexamples for the liquid crystal 
system (jl.ip . A very interesting open question is whether there are finite time singularities 
for the incompressible liquid crystal flows in two dimensions with finite energy and nontrivial 
velocity. We note that such nontrivial counterexamples (if they exist) must be non-radial 
since the incompressible constraint in the momentum equation of (11. ip destroys the spherical 
symmetry. 

2 We are grateful to Prof. Fanghua Lin for pointing this out to us. 
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By standard local theory (cf. [H] and [1'6\ I14j). there exists a local smooth solution to 
(jl.ip with the initial data f)1.2|) under the assumptions stated in Theorem II .li Moreover, by 
the weak-strong uniqueness (cf. [12] and [IS]), such a solution is also unique. For simplicity of 
presentation we will refrain from proving such local wellposedness results and only work with 
smooth local solutions. The main work then is to prove the a priori estimate (jl.5p . For this 
we will invoke the geometric angle assumption and prove that the corresponding harmonic 
energy is coercive. A crucial step is to establish the rigidity of approximate harmonic maps 
(cf. (|2.6p ) under the geometric angle condition. We state it as the following 

Theorem 1.5. Let e > 0, Cq > 0. There exists a positive constant 5q = 5o( e o, Co) such that 
the following holds: 

Ifd:R 2 ^ § 2 ; Vd G H 1 ^ 2 ) with \\Vd\\ L 2 < C and d 3 > e , then 

\\Vd\\i 4 <(l-5 )\\Ad\\ 2 L i. 

Consequently for such maps the associated harmonic energy is coercive, i.e. 

|| Ad + |Vd| 2 d|| 2 2 > | (|| Ad||| 2 + ||Vd||£ 4 ) . (1.8) 

Remark 1.6. Applying the method of [3] one can even obtain a stronger rigidity theorem. 
We thank Prof. Fanghua Lin for pointing this out. Our approach here is based on frequency 
localization combined with concentration-compactness. 

Note that our rigidity Theorem 11.51 is purely " static" and it has nothing to do with the 
velocity field of the fluid dynamics. As such it is stated for any map d : R 2 — >• § 2 under 
general assumptions. By using essentially the same arguments, we also have a similar result 
for the heat flows of harmonic maps: 

dt = Ad+ |Vd| 2 d, (1.9) 

where d is still a unit-vector on the sphere § 2 C R 3 (see [3] for a generalization of Eells- 
Sampson's theorem). We state the result as the following 

Theorem 1.7. Let Vd G L 2 (R 2 ) with \d \ = 1, d -e 3 G L 2 (R 2 ) and mf xeR 2 d 03 > 0. Then 
there exists a unique global smooth solution to the heat flows of harmonic maps (|1.9|) with 
the initial data d(0,x) = do(x). Moreover, (II. 7ft is satisfied. 

The proof of Theorem 11.71 will be omitted since it is essentially a repetition of the proof 
of Theorem 11.11 

We close this introduction by setting up some notations and conventions which will be 
used in this paper. 

Notations 

For any two quantities A and B, we use A < B (resp. A > B ) to denote the inequality 
A < CB (resp. A > CB) for a generic positive constant C. The dependence of C on other 
parameters or constants are usually clear from the context and we will often suppress this 
dependence. The value of C may change from line to line. For any function / : M 2 — > R, 
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we use HZ 1 1 lp or sometimes ||/|| p to denote the usual Lebesgue LP norm of a function for 

1 < p < oo. We use to denote the space-time norm 

ll/lk?^ = (/ (/ |/(i,z)Rz)")", 

with the usual modifications when p or r is infinity, or when the domain K x I 2 is replaced 
by a space-time slab. When p = r we abbreviate I? t U x by Lf x or L^. 

Occasionally we will need use the fractional Laplacian operator |V| S , s > 0, which is 
defined via Fourier transform as 

The homogeneous Sobolev norm H s for any s > is defined as ||/||u- s = |||V| S /||2 or more 
explicitly: 

11/11*3 = (/ \c\ 2s mm\ 2 dcy. 

We will need to use the Littlewood-Paley (LP) frequency projection operators. For sim- 
plicity we shall fix the notations on M. 2 , but it is straightforward to define everything in M. d 
for any d > 1. To fix the notation let <j) € C^°(M 2 ) and satisfy 

< 4> < 1, = 1 for |x| < 1, <j)(x) = for \x\ > 2. 

For two real positive numbers a < f3, define the frequency localized (LP) projection operator 
Pa<-<p by 

P a <.<pf = ~ 4>{a- l 0\Hf))- 

Here J- and J-~ l denote the Fourier transform and its inverse transform, respectively. Simi- 
larly, the operators P <a and P>^ are defined by 

P<ef = J r - 1 [ttP~ 1 $Mf)), 

and 

P >Q / = ^- 1 ([l-<A(a- 1 0W)). 
We recall the following Bernstein estimates: for any 1 < p < q < oo and dyadic N > 0, 

2 2 

\\P<Nf\\LUM*) ^ ^"'ll/lligCR 2 )- 

2 Rigidity 

In this section we prove our rigidity Theorem 11.51 Our proof is purely analytic, and uses 
in a quantitative way the geometric angle condition. The key ingredient of the proof is a 
frequency localization argument and a concentration-compactness procedure. 

We begin with the following simple lemma which locks the nontrivial L 2 weak limit of an 
L 2 bounded sequence of functions whose frequency is essentially localized to a (large) ring. 
The proof is stated for M. 2 but it naturally generalizes to M. d for any dimension d > 1. 
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Lemma 2.1 (Non-evacuation of L 2 -norm). Let Cq and a be two positive constants, and N 
be a dyadic integer. Suppose that f n : R 2 — > K is a sequence of functions such that 

\\fn\\L 2 < c o, Vn, 

and 

\\P± < ^ N fn\\L«' > a > 0, Vn. (2.1) 

jv ^ ^ 

Then one can find a subsequence (which we still label as f n ) and centers x n 6 M 2 , such that 

fn(x n ~ •) / weakly in L 2 , 

and 

\\fh*>l. 

Proof. By (|2.ip . we can find x n € M 2 such that 

ct 

\[Pj. < . <N fn){Xn)\ > 2> Vn - 

Or more precisely, for ip = J 7 ^ 1 ^(N^ 1 ^) — f(N£)^ £ <S(R 2 ), one has 

ct 

~4>(y)fn{Xn - y)dy\ > -, Vn. (2.2) 
z 

Since by assumption the sequence f n (x n — •) is uniformly bounded in L 2 , there exists a 
subsequence (we still denote it by f n {x n — -)) which converges weakly in 1? to some function /. 
Clearly, by taking the limit in (12. 2p and using the Cauchy- Schwartz inequality, one concludes 
that \\fWv > 1. □ 

Our next lemma allows us to remove the translational degrees of freedom in studying the 
rigidity property of approximate harmonic maps from M 2 to § 2 under the geometric angle 
condition. By removing the translational degrees of freedom (and quotienting out other 
possible non-compact group actions), we can restore the compactness in the same spirit as 
the usual concentration-compactness procedure. 

Lemma 2.2 (Removing translational invariance). Let Co and a be two positive constants. 
Suppose a sequence of maps d n : M 2 — > S 2 , n > 1 satisfy the following conditions: 

• \\Vd n \\ L 2 + HVdnlljji < C < oo, Vn, 

• ||Vd n || 2 4 > a > 0, Vn. 

Then, upon a subsequence in n if necessary, we can find a sequence of points x n 6 M. 2 such 
that 

\Vd n (x n — -)| 2 — /, weakly in L 2 

and \\f\\ L 2 > 1. 
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Proof. Let N be a dyadic integer which will be taken sufficiently large in the course of the 
proof. Our main idea is to localize the Vd n -sequence within the frequency window [1/N,N]. 
After that we apply Lemma 12.11 to find the nontrivial weak limit. 

We first deal with the high frequency piece. By frequency localization, we have 

P>iv((VcQ 2 ) = P >N (P >N/8 Vd n ■ P >N/8 Vd n ) + 2P >N {P< N/8 Vd n ■ P >N/8 Vd n ). 

By Sobolev embedding, we then have 

|[P>jv((V(i n ) 2 )|| L 2 < \\P >N / 8 Vd n \\ L 4\\Vd n \\ L 4 

< |||v|5p >Ar/8 vcy L 2 • c 

< N- l ' 2 Cl < a/10, (2.3) 

if we take N large enough. 

Similarly for the low frequency part, we use Bernstein's inequality to get 

ll^((Vci n ) 2 )||L 2 < AT 1 !!^) 2 )^! < a/10, (2.4) 

where again we need to take iV large enough. 

Now using (|2.3p and (|2.4p and the assumption || (Vd^) 2 1|_^2 > a, we obtain 

ll^<.<iv((Vdn) 2 )||L 2 > |- 

On the other hand, by Holder, one has 

<ll(Vdn) 2 ||| 1 ||||P^ < . <W ((V ( i n ) 2 )||l 00 
<l|^<.<iv((^n) 2 )|l!oo. 

Obviously 

ll^<. <JV ((Vdn) 2 )||r«- >1- 

We can then apply Lemma 12.11 to the sequence f n = |Vd n | 2 to get the result. □ 

Now we are ready to prove the rigidity Theorem 11.51 under the geometric angle condition 
d3 > eo- 

Proof of Theorem \1.5[ First of all, using a scaling argument d(x) —¥ d(Xx) with A = ||Ad||7 2 , 
we may assume 

\\Ad\\ L * = 1. 

It then suffices to show there exists So = 6o(eo, Co) > such that 

||Vd||| 4 <l-<5o. (2.5) 
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Assume (|2.5|) does not hold. Then we can find a sequence d n : M 2 -> § 2 such that 

ll^^nlli 2 = 1; II Vein || L 2 ^ ^0, d„3 > 6o, 

but 



Denote 



||Vd n ||x4 / 1, as n — > oo. 
5n = A(i n + jVd n | 2 d n . (2.6) 



It is easy to compute that 

HAW His = HAdnllia + ||Vd n ||^ 4 + 2 J \Vd n \ 2 Ad n ■ d n dx 

= l-\\Vd n \\\ A . (2.7) 

Hence 

IbnllL 2 \ as n — y oo. 

Applying Lemma 12.21 and performing a simple transform to eliminate the translation invari- 
ance if necessary, we conclude that 

|Vd n | 2 — 1 g(x), weakly in L 2 as n — > oo 

for some g(x) > and HffHia > 1. Consequently, by taking the limit n — > oo for the equation 
of d n s in (j2.6|) . we have 



Ad* 3 + g(x)e < 0, 

where d 3 is the weak limit of the sequence (upon subsequence if necessary) d n3 which satisfies 

||Ad*|| L2 < 1, \\Vdl\\ L 2<C . 

Now c?3 is superharmonic on R 2 and is bounded. Obviously d 3 must be a constant map. But 
since Ad^ + g(x)eo < 0, one concludes that g(x) must be identically 0. This contradicts to 
the fact ||<7||i,2 > 1. 

We have obtained the desired contradiction and therefore (|2.5p holds. Finally (ll.8p follows 
from (|2.5p . and the same computation (with d n replaced by d) as in (|2.6p and (|2.7p . 

This concludes the proof of Theorem 11.51 □ 

3 Proof of Theorem 11.11 

As was already mentioned in the introduction, we only need to establish the a priori estimate 
(jl.5p . By using the maximum principle (applied to the third component of d$), we have 

inf d 3 (t,x) > inf (fo3, Vt > 0. 



8 



One can apply the rigidity Theorem 11.51 to conclude that 

||Ad(*) + |Vd(t)| 2 d(t)||| 2 > §(||Ad(i)|| 2 L2 + ||Vd(i)||t) . V / > 0. (3. 1 ) 



where 5q > is the same constant as in Theorem 11,51 
Consequently, the basic energy law (|1.4p yields 



l\Nd(t)\\h + j£ (IIM*)lli» + \\Vd(t)\\h)dt (3.2) 



<~(IKII| 2 + ||Vd ||£ 2 ), Vt>0. 

This clearly settles (|1.5p and the first two terms in the estimate (jl.7p . It remains for us to 
prove the estimate 

roo 

/ \\V d{t)f L2 dt <oo. (3.3) 
J o 

By performing the basic energy estimate for the equation of d — e 3 and using the inter- 
polation inequalities, we have 

- 2 j t \\d-e,f L , + \\Vdf L2 

< J \Vd\ 2 d(d-e 3 )dx< ||Vd||2 4 ||d- e3 || i2 

< HVdlliallAdllylld-esll^ 
<^\\Vd\\l 2 +C\\Ad\\ 2 L2 \\d-e 3 \\ 2 L2 , 

where C > is an absolute constant. A Gronwall in time argument then gives 

\\d(t,-)-e 3 \\ 2 L2 + \\Vd\\ 2 L 2ds< \\d — e 3 \\\ 2 exp ( / || AdW^ds^j , Vi > 0. 

Noting (|3.2p . one concludes that 

ft 

||Vd||| 2 ds < 00, Vt>0. 



/ 

Jo 



This obviously implies 

This concludes the proof of Theorem ll.il Finally we point out that a further higher order 
energy estimate also implies that any H s norm of (it, Vd) is uniformly bounded in time for 
any s > 0. 
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